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Abstract 

We demonstrate that if the dark matter (DM) in the Universe contains multiple 
components, the possible interactions between the DM components may convert 
the heavier DM components into the lighter ones. It is then possible that the 
lightest DM component with an annihilation cross section significantly larger than 
that of the typical weakly interacting massive particle (WIMP) can obtain a relic 
density in agreement with the cosmological observations, due to an enhancement 
of number density from the DM conversion process at late time after the ther- 
mal decoupling, which may provide an alternative source of boost factor relevant 
to the positron and electron excesses reported by the recent DM indirect search 
experiments. 
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1 Introduction 



In the recent years, a number of experiments such as PAMELA pQ , ATIC [2] , Fermi-L AT 
[3] and HESS p[| etc. have reported excesses in the high energy spectrum of cosmic-ray 
positrons and electrons over the backgrounds estimated from the traditional astrophysics. 
Besides plausible astrophysical explanations [SHI], the dark matter (DM) annihilation or 
decay provides exciting alternative explanations from particle physics. 

If the DM particle is a thermal relic such as the weakly interacting massive particle 
(WIMP), the thermally averaged product of its annihilation cross section with the rela- 
tive velocity at the time of thermal freeze out is typically (av)p ~3x 10 _26 cm 3 s _1 . The 
positron or electron flux produced by the DM annihilation can be parametrized as 

$ e = BN e P ^, (1) 

where po — 0.3 GeV • cm -3 is the smooth local DM energy density estimated from 
astrophysics, N e is the averaged electron number produced per DM annihilation which 
depends on DM models and propagation parameters, and is the mass of the DM 
particle. The boost factor B is defined as B = (p/ p Q ) 2 {av) / (av) p with p the true 
local DM density and (av ) the DM annihilation cross section multiplied by the relative 
velocity and averaged over the DM velocity distribution today. Both the PAMELA and 
Fermi-LAT results indicate that a large boost factor is needed [HIE]. For a typical DM 
mass of ~1(1.6) TeV the required boost factor B is ~ 500(1000) for DM annihilating 
directly into p + p~ and p fixed to po [9]. 

A large boost factor may arise from the non-uniformity of the DM distribution in the 
DM halo. The N-body simulations show however that the local cumps of dark matter 
density are unlikely to contribute to a large enough p/p [10l[TT]. An other possibility 
of increasing the boost factor is that the DM annihilation cross section may be velocity- 
dependent which grows at low velocity. The DM annihilation cross section today may be 
much larger than that at the time of thermal freeze out, and thus is not constrained by 
the DM relic density. Some enhancement mechanisms have been proposed along this line, 
such as the Sommerfeld enhancement [T2H20] and the resonance enhancement [2X1 - 123"] etc. 

In some non-thermal DM scenarios, the number density of the DM particle can be 
enhanced by the out of equilibrium decay of some heavier unstable particles if the DM 
particle is among the decay products of the decaying particle j2U|25]. The decay of the 
unstable particle must take place at very late time. Otherwise the DM particles with the 
enhanced number density will annihilate into the Standard Model (SM) particles again, 
which washes out the effect of the enhancement. This requires that the decay width of 
the unstable particle must be extremely small, typically 10~ 17 GeV for the mass of the 
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decaying particle around TeV [23], which is much smaller than that of the typical weak 
interaction. 

In this work, we consider an alternative possibility for generating a boost factor, 
which does not require the velocity-dependent annihilation cross section or the decay of 
unstable particles. We show that in the scenarios of interacting Milt-component DM, 
the interactions among the DM components may convert the heavier DM components 
into the lighter ones, which is not sensitive to the details of the conversion interaction. 
If the interactions are strong enough and the DM components are nearly degenerate in 
mass, the conversion can enhance the number density of the lighter DM components at 
late time after the thermal decoupling. Eventually, the whole DM today in the Universe 
can be dominated by the lightest DM component with enhanced number density, which 
corresponds to a large boost factor. The scenarios of multi-component DM have been 
discuss previously in Refs. [26H35] . Note however that the models with simply mixed 
non-interacting multi-component DM cannot generate large boost factors. 

This paper is organized as follows: in section |2j we first discuss the thermal evolution 
of the DM number densities in generic multi-component DM models. We then give 
approximate analytic expressions as well as precise numerical calculations of the boost 
factor in a generic two-component DM model. In section [3j we consider a concrete model 
containing two fermionic DM particles with extra U(l) gauge interactions in the hidden 
sector. The conclusions are given in section HJ 

2 Thermal evolution of the interacting multi-component 
DM 

Let us consider a generic model in which the whole cold DM contains N components 
Xi (i — 1, ■ ■ ■ , N), with masses m 8 and internal degrees of freedom gi respectively. The 
DM components are labeled such that rrii < mj for i < j, thus xi is the lightest DM 
particle. We are interested in the case that Xi are nearly degenerate in mass, namely 
the relative mass differences between Xi an d Xi satisfy £j = (m.j — m 1 )/m 1 <C 1. In 
this case, we shall show that the interactions between the DM components lead to the 
DM conversion. The situation is analogous to the neutral meson mixing and neutrino 
oscillations in particle physics. They all occur at small mass differences. The thermal 
evolution of the DM number density normalized to the entropy density Yi = rii/s with 
respect to the rescaled temperature x = m%/T is govern by the following Boltzmann 
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equation 



dYjjx) A_ 
dx x 2 



.1 



(2) 



where A = xs/H(T) is a combination of x, the entropy density s and the Hubble param- 
eter H{T) as a function of temperature T. Y ieq ~ (gi/ s)[miT / (2ir)] 3 ^ 2 exp(— Six) is the 
equilibrium number density normalized to entropy density for non-relativistic particles. 
((Jiii) are the thermally averaged cross sections multiplied by the DM relative velocity 
for the process x%X% ~^ XX' with XX' standing for the light SM particles which are in 
thermal equilibrium, and (cr^v) are the ones for the DM conversion process x%X% ~~ ^ XjXj- 
The quantity 



Y ieq (x) ( g { \ ( m, 



3/2 

is the ratio between the two equilibrium number density functions, In writing down Eq. 
(J2J) we have assumed kinetic equilibrium. The first term in the r.h.s. of Eq.([2]) describes 
the change of number density of Xi due to the annihilation into the SM particles, while 
the second term describes the change due to the conversion to other DM particles. 

In the case that the cross section of the conversion process ((Tijv) is large enough, the 
DM particle Xi can be kept in thermal equilibrium with Xj f° r a l° n g time after both x% 
and Xj have decoupled from the thermal equilibrium with the SM particles. In this case, 
the number densities of Xi,j satisfy a simple relation 

Yiix*) ^ Yi e q(yX) f--(x\ (4) 

Yj\x) Yj eq (x) 

We emphasize that even when Xi is m equilibrium with Xj the ratio of the number 
density Yi(x)/Yj(x) can be quite different from unity and can vary with temperature. 
For instance, if gi 3> gj and < (e^ — €j) 1, from Eq. and (j4j) one obtains 
Yi(x) 3> Yj(x) at the early time when (ej — €j)x 1. However, at the late time when 
(ej — 6j)x ^> 1, one gets Yj,(x) Yj(x), which is due to the Boltzmann suppression factor 
exp[— (ej — €j)x] in the expression of r^. Thus the heavier particles can be gradually 
converted into lighter ones through this temperature-dependent equilibrium between x% 
and Xj- 

Since all the DM components Xi are stable, in general the co-annihilation process 
XiXj ~^ are n °t allowed as the crossing process Xi ~ > XjXX' corresponds to the de- 
cay of Xi- Furthermore, unlike the case of co- annihilation, x% an d Xj ma y n °t necessarily 
share the same quantum numbers. 

An interesting limit to consider is that the rates of DM conversion are large compared 
with that of the individual DM annihilation into the SM particles, i.e. ((Tijv) > (<Tiv). 
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In this limit, after both the DM particles have decoupled from the thermal equilibrium 
with the SM particles, which take place at a typical temperature x = Xd ec ~ 25, the 
strong interactions of conversion will maintain an equilibrium between Xi an d Xj f° r a 
long time until the rate of the conversion cannot compete with the expansion rate of the 
Universe. Making use of Eq. (TJJ, the evolution of the total density Y(x) = Y^h=i Y%{x) 
can be written as 

£ = -p<<V/»>0"-O. (5) 

where (a e ffv) is the effective thermally averaged product of DM annihilation cross section 
and the relative velocity which can be written as 

(a eff v) = — ^ ^ W}> ( 6 ) 

9eff 

where Wi = (<7ji>) / (<Jiv) is the annihilation cross section relative to that of the lightest 
one. The total equilibrium number density can be written as 

± f mi T\ i/2 



Y P „ = ^Y ieq (x) « g eff (^—) exp(-x), (7) 



eq — 

i=l 



with effective degrees of freedom g e ff = J2i9i(^- + £ i) 3 ^ 2 ex P( — Note that the con- 
version terms do not show up explicitly in Eq. ([5]). Through the conversion processes 
XiXi ~^ XjXj the slightly heavier components will be converted into the lighter ones, 
because the factor rij(x) is proportional to exp[— (m; — mj)/T] which suppresses the 
density of the heavier components at lower temperature. If the conversion cross section 
is large enough, most of the DM components will be converted into the lightest Xi before 
the interaction of conversion decouples, which may result in a large enhancement of the 
relic density of Xi an d leads to a large boost factor. 

As an example, let us consider a generic DM model with only two components. For 
relatively large conversion cross section u = (<J2\v) / (&iv) > 1, The effective total cross 
section is given by 

1 + w g 2 exp(— 2ex) 
[1 + gexp(— ex)f 

where w = w 2 , g = ^2/^1 and e = e 2 . Because of the x-dependence in (a e ffv), the 
thermal evolution of Y(x) differs significantly from that of the standard WIMP. In the 
case that X2 has large degrees of freedom but a small annihilation cross section, namely 
g ^> 1, w <C 1 and wg 2 1, the thermal evolution of the total density Y can be 
approximated by 

dY X 1 

dx x 2 [1 + g exp(-ex)] 2 x 



^effV) = h ; (CiV), (8) 



(Y 2 ~ Y e 2 q ) , (9) 
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the thermal evolution of the total number density can be roughly divided into four stages: 
i) At high temperature region where 3 < x <C Xd ec , both the DM components are 
in thermal equilibrium with the SM particles. Y{(x) must closely track Yi eq (x) which 
decrease exponentially as x increases. However, since g> 1 and e< 1, the number den- 
sity of X2 is much higher than that of Xi, i- e - Y 2 (x) ^> Yi(x). ii) When the temperature 
goes down and x is close to the decoupling point Xd ec , both the DM components start to 
decouple from the thermal equilibrium. In the region Xd ec ^ x <C 1/e, (a e ffv) is nearly 
a constant and (a e ffv) ~ (a\v)/{l + g) 2 <C {(T\v), the total density Y(x) behaves just 
like that of an ordinary WIMP which converges quickly to Y(x) ~ Xdec/i^i&iv ))• hi) 
As x continues growing, the suppression factor exp(— ex) in (a e ffv) becomes relevant. 
The value of (a e ffv) grows rapidly especially after x reaches the point ex ~ 0(1), which 
leads to the further reduction of Y(x). In this stage, although both xi,2 have decoupled 
from the thermal equilibrium with the SM particles. The strong conversion interaction 
X2X2 "H- X1X1 maintains an equilibrium between the two DM components. According 
to Eq. @, the relative number density Y 2 (x)/Yi(x) decreases with x increasing, which 
corresponds to the conversion from the heavier DM component into the lighter one. At 
the point x c = {1/e) In g one has Y 2 {x) ~ Y\(x). For the region x > Xd ec and x is not 
close to x c , because of Y eq (x) <C Y(x) and gexp(-ex) ^> 1, the Eq. d^J) can be analyt- 
ically integrated out, using the expression I(x) = f x~ 2 exp(x)dx = Ei(x) — exp(x)/x 
where Ei(rr) is the exponential integral function. The integral has an asymptotic form 
of I(x) ~ exp(x)/x 2 for x ^> 1. Thus Y(x) in this region can be approximated by 



iv) When x becomes very large ex ^> 0(1) , (a e ffv) quickly approaches (criv), and 
becomes independent of x again. The evolution of Y(x) in this region can be obtained 
by a simple integration as it was done in the stage ii). The solution of Y{x) shows a 
second decoupling. Finally when the conversion rate cannot compete with the expansion 
rate of the Universe at some point xp corresponding to sY 2 (a2iv) / H m 1, both Yi(x) 
and Y 2 (x) remain unchanged as relics. The whole DM can be dominated by xi if the 
conversion is efficient enough. 

By matching the analytic solutions of Y(x) in different regions near the points Xd ec 
and x c , and requiring that the final total relic density is equivalent to the observed 
^CDMh 2 ~ 0.11, we obtain the following approximate expression of the boost factor 



As expected, the enhancement essentially comes from the conversion of the degrees of 
freedom. Thus the maximum enhancement is g 2 . The two terms in the r.h.s of the above 




(10) 




(11) 
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equation correspond to the reduction of Y(x) during the late time conversion stages. For 
large enough g, the boost factor can be approximated by B « g 2 / (1 + eg 2 Xd ec / fog). In 
order to have a large boost factor, a small e <C In g/(g 2 Xdec) is also required. As shown 
in Eq. (TlTTl the boost factor is not sensitive to the exact values of the cross sections as 
long as the conditions w C 1 and u ^> 1 are satisfied. 

We numerically calculate the thermal evolution of Y^ix) and the boost factor without 
using approximations for a generic two-component DM model. The results for w = 10 -4 , 
u — 10 and e = 2 x 10~ 4 is shown in Fig. [TJ The value of (o^i?) is adjusted such that 
the final total DM relic abundance is always equal to the observed value QcDMh 2 . The 
mass of the light DM particle is set to mi = 1 TeV. For an illustration the ratio between 
the internal degrees of freedom is set to be large g = 60. From the figure, the four stages 
of the thermal evolution of Y(x) as well as the crossing point can be clearly seen. The 
crossing point at x = x c ~ 2 x 10~ 4 indicates the time when the number density of xi 
start to surpass that of \2 and eventually dominant the whole DM relic density. In this 
parameter set a large boost factor B « (civ ) / (av) p « 585 is obtained which is in a 
remarkable agreement with Eq. ( flTT) with error less than ~ 5%. For a comparison, in 
Fig. [TJwe also show the cases without conversions. 

In Fig. [2] (left), we show how the boost factor B varies with the mass difference e for 
different relative internal degrees of freedom g. In general, B becomes larger for smaller 
e and larger g. For e = 10~ 4 and g = 60, the boost factor can reach B ~ 10 3 . For a 
much smaller g = 20 and a larger e = 8 x 10~ 4 , the boost factor can still reach 0(100). 
The dependence of B on the cross sections u and w is shown in Fig. [2] (right). A small 
w and large u lead to the increasing of B. However, for very small w < 10~ 4 and vary 
large u > 100, the value of B becomes insensitive to the exact values of w and u, which 
is also in agreement with the approximate solution given in Eq. ( TTTT) . 

3 A simple model with DM conversion 

For models with multiple DM components, it is nature that there exists interactions 
among the DM components which may lead to the conversions among them. In this work 
we consider a simple interacting two-component DM model by adding to the standard 
model (SM) with two SM gauge singlet fermionic DM particles Xi,2- The particles 
Xi,2 are charged under a local U(l) symmetry which is broken spontaneously by the 
vacuum expectation value (VEV) of a scalar field <fi through the Higgs mechanism. The 
corresponding massive gauge boson is denoted by A which may cause the interaction 
X2X2 XiXi- The stability of Xi,2 is protected by two different global U{1) number 
symmetries. An SM gauge singlet pseudo-scalar r\ is introduced as a messenger field 
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Figure 1: Thermal evolution of the number densities Y\(x) (red solid) and Y2 (x) (blue 
solid) with respect to x. The solid (dashed) curves correspond to the case with (without) 
DM conversions. The green dotted curve corresponds to the sum of Y\ and Y 2 , for 
parameters g = 60, mi = ITeV, e = 2 x 10~ 4 , w = 10 -4 and u — 10 respectively. 

which couples to both the dark sector and the SM sector. In order to have the leptophilic 
nature of DM annihilation, we also introduce an SM SU(2)l triplet field A with the SM 
quantum number (1,3,1) and flavor contents A = (5 ++ ,5 + ,5°). The triplet carries the 
quantum number B — L=2 such that it can couple to the SM left-handed leptons 
through Yukawa interactions £ c l A£l, but cannot couple to quarks directly. The VEV of 
the triplet has to be very small around eV scale, which is required by the smallness of 
the neutrino masses. As a consequence, the couplings between one triplet and two SM 
gauge bosons such as S ±:iz W T W T , 5 ± W T Z° and 5°Z°Z° are strongly suppressed as they 
are all proportional to the VEV of the triplet, which makes it difficult for the triplet to 
decay even indirectly into quarks through SM gauge bosons [361137] . If rj has a stronger 
coupling to A than that to the SM Higgs boson H and <fi then the annihilation products 
of the dark matter particles Xi,2 will be mostly leptons. 

The full Lagrangian of the model can be written as C = Csm + £1 The new in- 
teractions in C\ which are relevant to the DM annihilation and conversion are given 
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Figure 2: Left) boost factor B as a function of the relative mass difference e for different 
relative degrees of freedom g=80 (solid), 60 (dashed), 40 (dotted) and 20 (dot-dashed) 
respectively, for w = 10~ 4 and u = 10 2 ; Right) boost factor as a function of the relative 
conversion cross section u. Four curves correspond to w = 10~ 5 (solid), 10~ 4 (dashed), 
5 x 10~ 4 (dotted) and 10~ 3 (dot-dashed) respectively, for parameters g = 60, ni\ = 1 
TeV and e = 1 x 10~ 4 respectively. 

by 

A D Xi(z# - Tm)xi + (D^)\D^) - m 2 ^4> 

+ ^d^r]d^r] - ^m 2 v r] 2 - yoCiijsVXi ~ ydl^h + h.c 

-(.VV + irf) [Tr(A f A) + k(H^H) + C(0V)] > («' = 1, 2) (12) 

with = d^ + igA^-fji and standing for the gauge coupling constant. Note that and 
77 do not directly couple to the SM fermions. After the spontaneous symmetry breaking 
in V ((/>), the scalar (J) obtains a nonzero VEV ((f)) = which generates the mass 

of the gauge boson vtla = Qav^- At the tree level, the three components of the triplet 
S ++ , 5 + and 5° are degenerate in mass, i.e. ms++ = ms+ = ms+ = m&. 

After the spontaneous symmetry breaking in the scalar sectors, the fields A, H and <f> 
obtain nonzero VEVs, which also generates a linear term in 77 through the last term of Eq. 
( TT2l) . The linear term in 77 in turn leads to a nonzero VEV of 77, i.e., (77) = v v ^ 0, which 
will give corrections to the masses of Xi an d may enlarge the mass difference between xi 
and X2- This problem can be avoided by using the above mentioned assumption that rj 
has a much stronger coupling to A than that to H and 0, which requires that k, C <C 1. 
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The VEV of r\ is proportional to the ratio between the linear and quadratic terms in 77, 
and can be estimated as v v ~ —h(kv 2 h + (v 2 ) / (2(m 2 + nv 2 H + (vl)). Since the VEV of 
the triplet A is extremely small and Vh ~ O(10 2 )GeV, if /1, m v , and are all around 
TeV scale, for k < C(10" 2 ) and C < C(10" 4 ) the VEV of 77 is v v < C(10~ 4 )TeV which 
is small enough to avoid breaking the degeneracy in the masses of Xi,2- 

We assume that \2 has large internal degrees of freedom relative to that of Xi, i- e -; 
92 3> 9i, which can be realized if X2 belongs to a multiplet of the product of some global 
nonabelian groups. For instance #2 = 4#2 with c/2 =16, 8, and 4 if it belongs to the 
spinor representation of a single group of SO (8), SO(Q) and SO (4) respectively. When 
X2 belongs to a representation of the product of these groups, its internal degrees of 
freedom can be very large. 

At the early time when the temperature of the Universe is high enough, the triplet 
A can be kept in thermal equilibrium with SM particles through the SM gauge interac- 
tions. The DM particles x% can reach thermal equilibrium by annihilating into the triplet 
through the intermediate particle 77. The annihilation X2X2 — > V* ~ >* 5 ±± 5 =FT , 5 ± 5 T , 5°5°* 
is an s-wave process which is dominant contribution . The cross section before averaging 
over the relative velocity v is given by 



where Nf = 3 is the number of final states, m v is the mass of 77 and s is the square 
of the total energy in the center of mass frame. For s-wave annihilation we use the 
approximation that the thermally averaged cross section is the same as the one before 
the average, i.e., (av) ~ av. From the above equation the ratio of the two annihilation 
cross sections is w = {1/2 1 Hi) 2 {gi 1 92)- It is easy to get a very small w provided that 
V2 <S yi and gi g 2 . In order to have a large enough (aiv) ^> (ctv)f the product of 
the coupling constants y\\i must be large enough, or the squared mass of 77 is close to s. 
The cross section of the conversion process X2X2 — t A* — > X1X1 is given by 



s is close to 4m 2 . at the vary late time of the thermal evolution. However, the cross 
section be greatly enhanced if rriA is close to a resonance when the relation s ~ m\ is 
satisfied. In the numerical calculations, we find that for the following selected parameters: 
mi = ITeV, e = 1 x 10~ 4 , g x = 1, g 2 = 60, m A = 500 GeV, m v = 1.5 TeV, m A = 2.02 
TeV, yi — 3, y 2 — 0.07, fi/rrii = 3, and g& = 2.5, the following ratio of the cross section 
can be obtained 




(13) 




(14) 



The cross section is suppress by g\ / ^2 and also the phase space factor a/1 — Amf/s when 



w ~ 1 x 10 5 , u ~ 0.5, and (<J\v) / '{crv) f — 500. 
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In this parameter set the relative mass difference between rriA and 2r«2 is around 1%. 
From Fig. [21 one can see that the corresponding boost factor is B ~ 500, which is large 
enough to account for the PAMELA data for the dark matter mass around TeV. 

4 Discussions and Conclusions 

The mechanism proposed here does not require velocity- dependent annihilation cross 
sections which is essential to the Sommerfeld enhancement. There exists stringent con- 
straints from astrophysical observations if the DM annihilation cross section scales with 
velocity as 1/v or 1/v 2 and saturates at very low velocity. Those constraints involves 
the bound on the /i-type distortion of CMB spectrum [38H1U] and the bounds on diffuse 
gamma-rays from the cold structures which have lower velocity dispersion than that in 
the solar neighborhood in which v ~ 10~ 3 . For instance, in the subhalos the average 
velocity can be as low as v ~ 10~ 5 [41], and the DM velocity in the protohalos can 
be even lower v ~ 10 -8 [12] . If the enhancement is insensitive to the velocity, those 
apstrophysical bounds can be relaxed significantly. Furthermore, unlike the Sommerfeld 
enhancement, no attractive long-range force between the DM particles is involved. The 
existence of such a long-range force can change the halo shape and is constrained by 
observations [43H45] . The boost factor from DM conversion is free from this type of 
constraint as well. 

In summary, We have considered an alternative mechanism for obtaining boost fac- 
tors from DM conversions which does not require the velocity-dependent annihilation 
cross section or the decay of unstable particles. We have shown that if the whole DM 
is composed of multiple components, the relic density of each DM component may not 
necessarily be inversely proportional to its own annihilation cross section. We demon- 
strate the possibility that the number density of the lightest DM component with an 
annihilation cross section much larger than (av) p can get enhanced in late time through 
DM conversation processes, and finally dominates the whole relic abundance, which cor- 
responds to a boost factor needed to explain the excesses in cosmic-ray positron and 
electrons reported by the recent experiments. 
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